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Abstract 

We study the inflation scenarios, in the framework of superstring theory, where the inflaton 
is an axion producing the adiabatic curvature perturbations while there exists another light 
axion producing the isocurvature perturbations. We discuss how the non-trivial couplings 
among string axions can generically arise, and calculate the consequent cross-correlations 
between the adiabatic and isocurvature modes through concrete examples. Based on the 
Planck analysis on the generally correlated isocurvature perturbations, we show that there is 
a preference for the existence of the correlated isocurvature modes for the axion monodromy 
inflation while the natural inflation disfavors such isocurvature modes. 
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1 Introduction 


The recent Cosmic Microwave Background (CMB) data is still consistent with the simple ACDM 
model with a nearly scale-invariant, adiabatic and Gaussian power spectrum which can well 
be represented by the single-held slow-roll inhation models HHa. The forthcoming cosmolog¬ 
ical data with even better precision however could reveal the potential deviations from such 
pure adiabatic perturbations, and it would be worth exploring the possibilities for the non- 
adiabatic perturbations in existence of correlations among the adiabatic and non-adiabatic 
modes along with their indications for the early Universe dynamics in the fundamental physics 
models. In this paper, we study the mixture of the adiabatic mode and cold dark matter (ax- 
ion) isocurvature mode taking account of their possible cross-correlations [1HI2] through the 
concrete models based on superstring theory. The recent Planck analysis studied the generally 
correlated isocurvature perturbations. The robust parameter estimation, without signihcantly 
affecting the bounds on the conventional (ACDM) cosmological parameters even with the in¬ 
clusion of isocurvature modes, was not previously realized due to strong degeneracies among 
the parameters involving the isocurvature perturbations which the previous CMB data sets 
suffered from and the Planck data could greatly reduce [3]. With such a precise cosmological 
parameter estimation including the correlated isocurvature perturbations at hand, it would be 
intriguing to explore the indications of the generally correlated isocurvature perturbations for 
the early Universe phenomena, and we in this paper aim to study the presumably ubiquitous 
light degrees of freedom in the early Universe through their isocurvature fluctuations. 

In superstring theory, the higher-dimensional form fields predict a number of light degrees 
of freedom represented by the axions in addition to the QCD axion [13]. These axions are 
associated with the internal cycles of the extra-dimensional manifold. While, at the pertur¬ 
bative level, the axion potential is protected by the gauge symmetry in string theory, the 
non-perturbative effects can break the continuous gauge symmetry leading to the discrete one 
and generate the axion potential. It is thus expected that the axion potential is well controlled 
by the residual discrete symmetry and the mass scale of axions depends on the non-perturbative 
effects [HHTB]. We, in this paper, focus on the single-held axion inhation models in coexistence 
of an isocurvature perturbation due to another light axion. 

Although the axion inhation is often considered as the single-held inhation, the axion po¬ 
tential, in general, has the axionic mixing due to the moduli-mixing gauge kinetic function. 
As concrete examples, we discuss the natural inhation ini in Sec. [2] and axion monodromy 
inhation [TSlfTO] in Sec. [3] where the cross-correlated isocurvature perturbations can arise due to 
such axionic mixings. In these illustrative examples, the adiabatic curvature perturbations are 
dominantly sourced by the axion-inhaton whereas the isocurvature perturbations originate from 
the huctuations of the light axion (diherent from the heavy axion inducing the inhation). The 
mixing of the string axions arises from the non-perturbative ehects in the sinusoidal form, and 
the consequent cross-correlations between the adiabatic and isocurvature modes are studied. 
We conclude our discussions in Sec. 01 
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2 Natural inflation with sinusoidal correction 


The natural inflation is among the simplest axion inflation models UTl which can be constructed 
in the held theory as well as superstring theory. Although, at the perturbative level, the axion 
potential is not generated due to the gauge symmetry in string theory, the non-perturbative 
effects in a hidden gauge sector can generate the axion potential terms. Especially, when 
the gauginos (A) of the hidden gauge group condensate at a certain energy scale [21], the 
superpotential can be generated in the four-dimensional (4D) J\f = 1 supersymmetry, 

W ~ (AA) ~ (1) 

where A = 0(1) and a = 247r^/6o with Bq being the one-loop beta-function coefficiently] We 
consider the scenarios where the size of gauge coupling is determined by the real part of modulus 
field, T, which is typical for heterotic string theory, type I string theory and type II string theory 
with D-branes along the single cycle (see for reviews, e.g.. Refs. [m[T3] b By fixing the real part 
of moduli, for instance through another non-perturbative effect, we can obtain the effective 
inflaton potential for the imaginary part of modulus (axion), 

Rinf = At - cos|^ , (2) 

with (j) and / being the axion-inflaton and its decay constant. The conventional natural inflation 
model, in view of the recent Planck data, requires the trans-Planckian axion decay constant, 
/ > 5, even though its construction requires some care because the fundamental axion decay 
constant obtained after the dimensional reduction is typically much smaller than the Planck 
scale [22]. The gauge couplings in the visible and hidden sectors, in general, depend on the 
linear combination of moduli fields through the gauge threshold correction and non-trivial brane 
configuration. For example, in type II string theory, the D-branes wrap the internal cycle of 
six extra-dimensional manifold and then the volume of this internal cycle is determined by the 
linear combination of moduli fields Tj where the number of moduli Tj is determined by the 
topology of the extra-dimensional manifold [HllTS]. Thus, the gauge coupling on Dp-branes is 
represented by the linear combination of them, 

(ciT-) = 4 (3) 

where c* are constant. Furthermore, if we consider the one-loop corrections for the gauge 
coupling, the superpotential also depends on the linear combination of moduli fields T and T', 

W = (4) 

where A = 0(1), a = 247r^/6o and d = 247r^/6o x b/ASn with b being the one-loop beta-function 
coefficient determined by massive modes [23] . The axion decay constant for the modulus T' here 
can be enhanced by the one-loop effect [2U[25]. Indeed, there are several scenarios to enhance 

^Here and in what follows, we employ the reduced Planck units, Mpi = 2.4 x 10^® GeV = 1. 
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the axion decay constant based on the moduli-mixing in the gauge kinetic function such as the 
alignment mechanism [26], N-flation m, kinetic mixing [28], the threshold correction [2HI25] 
and the flux-induced enhancements |29j . 

Since there exist, in general, ubiquitous axion fields in string theory, one can also expect 
that there are moduli-dependent correction terms in the potential. 


= + (5) 

where 0 and y represent an axion-inflaton and another light axion held. For the notational 
brevity, we in the following dehne the parameters 
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( 6 ) 


so that the total potential can be written as 


V = A^(l — coscr) -I- A 2 (1 — cos {^jJ -f 6)). 


( 7 ) 


We hereafter focus on the scenarios where the adiabatic perturbations are dominantly sourced 
by the axion-inhaton huctuation 6(j) and the additional axion huctuation 6x leads to the isocur¬ 
vature perturbations. 

Before discussing the cosmological perturbations and their indication for the model discrim¬ 
ination, we show an allowed parameter region for the spectral tilt of the adiabatic perturbations 
and the tensor-to-scalar ratio in both pure adiabatic (ADI) model and generally-correlated ADI 
-|- cold dark matter isocurvature (GDI) model. Fig. [1] shows that, for the natural inhation with 
isocurvature perturbations, the inclusion of a cross-correlated isocurvature mode tightens the 
constraints on the axion decay constant, 5 < / < 10, and the e-folding number, 60 < N. On 
the other hand, the axion monodromy inflation, to be discussed in the next section, except 
for the quadratic one can be better fitted by the Planck data by including the cross-correlated 
isocurvature mode. 

The degeneracies among the parameters involving the correlated isocurvature perturbations 
result in the shift in the best-£t parameters compared with those in the pure adiabatic model, 
even though the strong degeneracies such as that between the isocurvature perturbation ampli¬ 
tude and adiabatic perturbation spectral index which WMAP data had greatly suffered from 
reduced significantly in Planck TT + polarization data inEi. The constraints on r however 
turn out not to be significantly affected by the inclusion of the cross-correlated isocurvature 
modes partly because the Planck data including the polarization already gives sufficiently tight 
constraints on the isocurvature and tensor modes [5]. In the rest of the paper, we for simplicity 
do not consider the significant tensor contribution and we adopt the Planck likelihood analysis 
results without including r in the following discussions. 

We now discuss the cosmic perturbations for the axion fields, starting with the brief dis¬ 
cussions for the conventional curvature and isocurvature perturbations to set up our notations 
followed by the exploration on their cross-correlations along with their indication for the string 
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N=60, Natural Inflation 

- N=50, Natural Inflation 
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Figure 1: 68% and 95% confidence level constraints on the adiabatic spectral index and tensor 
to scalar ratio from Planck [3]. The hlled contours are for generally-correlated adiabatic and 
CDM (axion) isocurvature modes. The unhlled dashed contours are for the pure adiabatic 
model without the isocurvature perturbations. 


inflation model building [20]. The curvature and isocurvature perturbations in our scenarios 
are 


n = 



( 8 ) 


X 


do 


(9) 


with fla and being the axion and matter densities with respect to the critical density. The 
factor Qa/^m appears here because we are interested in the isocurvature perturbations between 
the radiation and the non-relativistic matter, and the non-adiabatic fluctuations arise solely 
from an axion which contributes to the total matter density with the fraction Qa/^m- The dot 
denotes the time derivative and the subscript o represents the background held values during 
the inhation and, in the following, we omit 6 representing the huctuations for the notational 
brevity when it is clear from the context. 
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The corresponding power spectra are given by [30 
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6 V/ 


( 10 ) 


( 11 ) 


The cross-correlation between the curvature and isocurvature perturbations can be ob¬ 
tained using the in-in formalism 1311. and the relevant interaction term in the Hamiltonian 
at the quadratic order, A 2 cos(' 0 o + do)S'ijjS9, leads to the isocurvature cross-correlation power 
spectrum [201I3T] . 
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cos(0o + 


( 12 ) 


where the numerical integral of the Hankel function gives a factor ~ —0.45 and an order of the 
Hankel function is taken as u^(^g) = ^9/4 — We plot the following cross correlation 

parameter 




Vc 

y/V-jiPx 


~ 4.2 




cos(0o + 9q) ~ 4.2 




cos (00 + 9o) 


0g 

967r2 ’ 


(13) 


by varying the axion decay constant / from 1 to 20 and the e-folding number N from 50 to 60 
in Fig. 121 Here, the prefactor of 0 q is set to be of order 0.001 and the power spectrum of the 

^This integral can be evacuated, at the leading order, at an arbitrary value of x as long as < ^ < 1 

with ^ being the typical size of the slow-roll parameter [32l|33] . 
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adiabatic curvature perturbations 


'P'R. — 




(14) 


is fixed to be As — H'^/{87i^e) ~ 2.2 x 10“® with e ~ 2/0 q being the slow-roll parameter 
at the pivot scale /c* = 0.05 Mpc“^. This figure also shows the Planck likelihood contours 
including the polarization data which greatly improve the constraints on the isocurvature per¬ 
turbations compared with the WMAP results [HIS] . The high-/ (/ > 30) T E, EE data turn out 
to drive the isocurvature cross-correlation towards a smaller value and disfavor the negative 
cross-correlations which would be allowed otherwise with the high-/ TT data [3]. We can find 
that the coefficient c in (3c = c4>l has to be of order less than 10“^ to be within 2 sigma and 
the axion decay constant / is constrained to the range between 5 and 10. The cross correlation 
parameter (3c is constrained to be —0.1 ^ (3c ^ 0.3, or, in terms of the parameters in the 
sinusoidal correction term (Eq. ([5])), to be within 


-0.1 < 4.2 


9i92 




2.2 X 10-9 


cos('0o + Oq) 


OOtt^ 


< 0.3. 


(15) 


Moreover, the following conditions are taken into account to justify our calculations: 


The adiabatic perturbations come from Ijnf and not from Idnt, that is, Idnf ^ Idnt- 


avinf 

d(j> 


> 


dVint 

d<j> 


• The infiaton dynamics is dominated by 0, i.e., 

• The quantum fluctuations of axions are not over-damped during the inflation, m^, <C H^. 

• The standard slow-roll conditions, e <C 1 and \r]\ <C 1. 

Some of the above conditions may be redundant depending on the parameter range of interest. 
In the light of these conditions, the cross-correlation parameter is bounded above by 


(3c - 4.2 
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< 2.1 
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(16) 


where |A2Cos('0o + 6'o)| |ldnf| is applied. For gi,g 2 < C>(1), the constraint of Eq. (fT6|) for (3c is 
automatically satisfied if Eq. (1T5|) is satisfied. The illustrative values of isocurvature parameters 
are listed in Tab. [Uby setting the typical values for the parameters in the scalar potential in Eq. 

Note that, although the axion decay constants are typically of order the grand unification 
scale (10^9 GeV) [22] and hence one may expect gi ~ 5 ^ 2 , the hierarchical values gi g 2 
{91 ^ 92 ) can well be realized for the axion y (0) by the non-perturbative effects through the 
(gauge) threshold correction (Eq. (|1|)). 

Next, we estimate the fraction of isocurvature perturbations 
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where the power spectrum of the adiabatic perturbation is fixed as in Eq. flT^ . whereas the 
power spectrum of the isocurvature perturbation is given by 



Then, the fraction of isocurvature perturbations 
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(18) 


(19) 


can give a sizable contribution to the cosmological observables as illustrated in Fig. [3]where the 
prefactor of 0^^ in Eq. flT^ is set to 1 and 10 for a varying /. A larger prefactor is preferred 
for a larger isocurvature contributions. The Planck bounds the uncorrelated axion isocurvature 
mode to (3iso ^ 0.038, whereas the inclusion of isocurvature cross-correlation results in the 
constraint 0.034 < (3iso ^ 0.28 at the 95% confidence level [3]. Tab. [T] summarizes the typical 
numerical values of parameters in the scalar potential (Eq. ([7])) which can realize a sizable 
fraction of isocurvature perturbations. 
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0.964 

10 

55 

10 -^ 

1 
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0.07 

0.964 

10 

55 

1 

10-2 

1 X 10-5 

0.02 

1/2 

2 

1 X 10-5 

7 X 10-5 

0.964 


Table 1: The typical numerical values of parameters, the e-folding number {N), the spectral 
index (n^), the fraction of isocurvature perturbation {f3iso) and the cross-correlation parameter 
{Pc) for the natural inflation with sinusoidal correction. 


Before concluding this section to move on to the discussion on the monodromy inflation, let 
us comment on the microscopic description about the correction term in the inflaton potential 
given by Eq. ([5]). Such a potential can be derived from the following Kahler and superpotential, 

K =-2 ln(Ti + fi) - ln(T 2 + fs), 

IT = Wo + (20) 

where Wq is the flux-induced constant term induced by the Gukov-Vafa-Witten superpotential 
hTflux = f G A Q, where G is the linear combination of Ramond-Ramond and Neveu-Schwarz 
three-form fluxes and G is the period vector in the framework of type II superstring theory. 
The second and third terms in the superpotential fl20|) denote the non-perturbative effects, such 
as the gaugino condensation terms, D-brane instanton and world-sheet instanton effects. 

Let us dehne the moduli as 


Ti — -l- ioi, 

T 2 = G + (21) 
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Figure 2\ VdVV rPi and the adiabatic spectral index Ug for the natural inflation with sinusoidal 
correction (68 % and 95 % CL contours are from Planck [3]). Vd\/Vn^ = c x 0 q for 
c = 10“^, 5 X 10“^ are shown for varying iV and / (the labels are in units of the reduced Planck 
mass). The anti-correlation cases (for c = —10“^, —5 x 10“^) are also shown with the dotted 
curves. 
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C ; 68% CL 
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-■ N=50, P|/PFi=10/(t)^ 
f=20 


- f=10 


- f=5 

- N=60, P|/Pp|=1/^2 

- ' N=50, P|/Pr=1/^^ 

- f=20 


- f=7 

- f=5 


Figure 3: do = Vxj iVn. + Vx) and Ug for the natural inflation with sinusoidal correction (68% 
and 95% CL contours are from Planck [3]). (/Siso, ng) are shown for Vx/Vn = c/d with c = 1,10 
for varying / and N. 


and assume that all the real parts of moduli Ti ^2 are stabilized at their minima and sufficiently 
















heavier than the remaining imaginary parts of Ti^ 2 - Then, from the four-dimensional scalar 
potential based on 4D M = 1 supergravity, 

V = e^{K^^DiWDjW -3\W\‘^), (22) 


where is the inverse of the Kahler metric Kjj = d‘^K/d^^, DjW = Wi + KjW, with 
Wj = dW/d^^ and Kj = dK/d^^, for = T^,T^. We can obtain the axion potential for 
ai ,2 by further assuming that some uplifting sector lifts up the scalar potential from the AdS 
vacuum to the dS one with a very small vacuum energy. 


1/ ~ A Aicos -h A 2 COS , (23) 

where A ~ —Ai — A 2 — —Ai is a constant and Ai ^2 depend on the vacuum expectation values 
of Reri^ 2 - The helds 0 and y are the canonically normalized axions. The kinetic terms of 01^2 
are extracted from the second derivatives of the Kahler potential with respect to the moduli, 

Kjjd^^d^^ = Kt,t, dTi + Kt,t, dT2df2. (24) 


As a result, the axion decay constants /i, gi ^2 of the canonically normalized axions 0 and y are 
given by 
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3 Axion monodromy inflation with sinusoidal correction 

We now discuss the axion monodromy inflation which offers another popular axion inflation 
scenario in string theory. Axion monodromy inflation is a successful large-held inhation in 
which the inhaton can move around its conhguration place on many cycles, and the held range 
of inhaton can be thus much larger than its fundamental period determined by the axion decay 
constant. 

The scalar potential for the axion monodromy inhation is represented by 

c = -1(84,? - Mfr, ( 26 ) 

where 0 is the axion originating from the higher-dimensional form helds, represents the 
energy scale and p is the fractional number which depends on the model in string theory |16] . 
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Let us consider the spacetime filling D5-brane in type IIB string theory [18]. The D5-brane 
wraps a certain internal two-cycle S 2 in the 6D compact space in addition to the 4D spacetime 
and its Dirac-Born-Inf eld action is given by 

Sd5 = /r, \5-/ ,N3 [ d®Cr^—det(Gafe -|- Bab), (27) 


where Qs is the string coupling, a' is the regge-slope, Gab, a,b = 0, 1, 2, 3,4, 5 is the pullback of 
the metric of the target space, Bab is the Kalb-Ramond field whose extra-dimensional component 
corresponds to the axion b = B 2 where B 2 is the Kalb-Ramond two-form. We here do not 
consider the magnetic flux background. 

After carrying out the dimensional reduction, the axion potential can be extracted as 


Reff ^ 


T 

{2'KYgs{a'Y 


Vl* + IP, 


(28) 


where T and / are some warp factors and the volume of two-cycle S 2 in string units (a' = 1). 
For a large field value of the infiaton b ':$> the potential reduces to a linear type. 


Rff ~ /o \5 / /\2 

{27rfgs{a'y 

Then, the relevant Lagrangian of the infiaton is given by 

C = - i44>, 


(29) 


(30) 


where jaf = f(^ 2 -K)^g {a')^ with / being the decay constant of the axion 0 = 6. Furthermore, for 
the D4-brane in a nilmanifold (twisted torus) on type IIA string theory, the axion potential has 
the form of Eq. fl26|) with p = 2/3 [19]. When we consider the seven-branes [M] or a four-form 
field strength [35], the axion monodromy inflation is that with p = 2. The other types of axion 
monodromy inflation with p = 4/3,3 can also be constructed by a coupling between NS-NS 
two-form and the Ramond-Ramond field strength [3B] . 

As mentioned for the natural inflation, the axion, in general, can receive the non-perturbative 
effects associated with the gaugino condensation, D-brane instanton and world-sheet instanton, 
and the scalar potential receives the moduli-dependent correction including the mixing with 
another light axion y. 


V = pi\ ^(jf + yU^cos (— + —^ , (31) 

\9i 92 J 

where g 2 denotes the decay constant of y. We here assume that the moduli except for the 
relevant axions under our discussion are fixed at their minimum and decoupled from our setup. 
For the notational brevity, we in the following define the parameters 
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so that the total potential can be written as 

V = + /I 2 cos (ip + 0). (33) 

Analogously to the natural inflation discussed in the last section, the curvature and isocur¬ 
vature perturbations in our monodromy inflation scenario read 


with 


-Pn = Ask^^-\ 
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(34) 


(35) 


replacing /ii by As through the CMB normalization, ^ where p = 2 

corresponds to the natural inflation case. The cross-correlation power spectrum then becomes 


Pc ~ 2.1 
= 2.1 


9192 
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We plot the cross-correlation parameter 
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(37) 


by varying the index p and the e-folding number N in Fig. 01 where the prefactor of (fi^/p^ 
is set to ±0.005 for concreteness, /^c increases for a larger p because the initial held value of 
axion-inflaton increases for a larger p. By considering the consistency conditions as spelled out 
below Eq. ffT5|) in the last section for the validity of our calculations, the cross-correlation is 
bounded above by 

(l^)^Fos(<f„ + »o) 

<2.l(^V (38) 

Tab. [2] lists some illustrative values for the isocurvature perturbation parameters by setting the 
typical values for parameter sets in the scalar potential (l3l|l . 

We next estimate the fraction of isocurvature perturbations 
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where the power spectrum of adiabatic perturbations is fixed to be Ag = 2.2 x 10 ® at the 
pivot scale /c* = 0.05 Mpc“^, whereas the fraction of Vn and Vx is given by 



Fig. Oplots /Siso with p = 2/3,1,4/3, 2 as a function 0o by setting the parameters as 

Vx _ fpV 

Vn '''' UoJ ’ 


(40) 


(41) 


with c = 1 , 10 , 50 , and we can find, as expected, the isocurvature contribution increases for 
a larger p. The Planck data hence favors the sizable generally correlated isocurvature per¬ 
turbations for the axion monodromy inflation with sinusoidal correction. Tab. [ 2 ] exemplifies 
the parameters which can realize the sizable fraction of isocurvature perturbations. Figs. [H 0] 
and 0] hence demonstrates that, for the axion monodromy inflation with p = 1 , 2/3 including 
the sinusoidal correction, there is a preference for the existence of cross-correlated isocurvature 
modes in the currently available CMB data. 
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Table 2: The typical numerical values for the axion monodromy inflation with sinusoidal cor¬ 
rection. 
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Figure 4: Vc / \/VtiPx and for the axion monodromy inflation (1/inf = with sinusoidal 

correction (68 % and 95 % CL contours are from Planck [3]). Vd\/VnPc = c x for 

c = 0.005 are shown for varying the e-folding number N. The anti-correlation cases (for 
c = —0.005) are also shown with the dashed lines. 
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Figure 5: /^iso = Vx/{V-jz + Vx) and n* for the axion monodromy inflation (Idnf = with 

sinusoidal correction (68 % and 95 % CL contours from Planck [3]). (Aso,ris) are shown for 
'Px/'P'R. = c X p^/(pl with c = 1, 5, 50. 
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4 Conclusion 


The natural inflation and axion monodromy inflation with moduli-dependent sinusoidal cor¬ 
rections can generically appear in the scalar potential through the non-perturbative effects. 
We demonstrated that probing the precise nature of isocurvature fluctuations can help us un¬ 
derstand the nature of fundamental physics using these popular inflation models as concrete 
examples. In this paper, we focused on the scenarios where the heavy axion induces the adia¬ 
batic perturbations while another light axion sources the isocurvature perturbations with their 
cross-correlations taken into account. 

Sec. [2] demonstrated that the cross-correlated isocurvature mode gives an even tighter con¬ 
straints on the decay constant of axion-inflaton and the e-folding number for natural inflation. 
While Sec. |2] showed that the cross-correlated isocurvature perturbations are not favored by 
Planck for the natural inflation, Sec. [3] showed that there is a preference for the existence of 
cross-correlated isocurvature mode for the axion monodromy inflation. 

We also mention that, when, in contrast to our setup, the sinusoidal corrections are not 
suppressed enough, the scalar power spectrum could posses the modulating behaviour for the 
natural inflation [2^H0] and the axion monodromy inflation [38l[39lHTlll2] . Such an additional 
feature in the cosmological observables would also be of great interest to discriminate among 
the possible inflation models. 

We illustrated our findings through the simple models where a single axion is added besides 
the axion-inflaton, but, in general, there could appear multiple light axions Xi addition to 
the axion-inflaton 0. They can then have the mixing terms 



and the analysis analogous to what has been done in this paper can be preformed for such 
multiple axion cases too. 

Even though Planck, in particular the addition of the polarization data, can significantly 
tighten the constraints on the isocurvature modes, we should be careful about the potential 
systematics in the current Planck data. For instance, it was pointed out that the apparent low 
power in rE spectrum could result in the preference for the positive cross-correlation, and this 
could cause an over-constraint on the isocurvature component if such a power spectrum feature 
is due to the unidentified systematics |13]. While the current polarization data at hand are 
not yet robust to the systematis on the large scales, the forthcoming polarization data with a 
better handle on the systematics would certainly be able to probe the nature of isocurvature 
perturbations more precisely and consequently explore the properties of ubiquitous light degrees 
of freedom in the early Universe. 
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